Quantum Zeno dynamics of a field in a cavity 
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We analyze the quantum Zeno dynamics that takes place when a field stored in a cavity undergoes 
frequent interactions with atoms. We show that repeated measurements or unitary operations 
performed on the atoms probing the field state confine the evolution to tailored subspaces of the 
total Hilbert space. This confinement leads to non-trivial field evolutions and to the generation of 
interesting non-classical states, including mesoscopic field state superpositions. We elucidate the 
main features of the quantum Zeno mechanism in the context of a state-of-the-art cavity quantum 
electrodynamics experiment. A plethora of effects is investigated, from state manipulations by phase 
space tweezers to nearly arbitrary state synthesis. We analyze in details the practical implementation 
of this dynamics and assess its robustness by numerical simulations including realistic experimental 
imperfections. We comment on the various perspectives opened by this proposal. 

PACS numbers: 03.65.Xp, 42.50.Dv, 42.50.Pq 



I. INTRODUCTION 

The evolution of a quantum mechanical system can be 
significantly slowed down by a series of frequent measure- 
ments pp. This effect, named after the Eleatic philoso- 
pher Zeno [2] , has attracted widespread attention during 
the last 20 years, since Cook proposed to test it on oscil- 
lating (two- level) systems [3]. This was a simplified ver- 
sion of the seminal idea by Misra and Sudarshan [2 , who 
had in mind genuinely unstable systems, but it had the 
important quality of making the quantum Zeno 'paradox' 
(as it was originally considered) amenable to experimen- 
tal test. 

The quantum Zeno effect (QZE) has been successfully 
demonstrated in many experiments on various physical 
systems, such as r.f. transitions between ionic hyperfine 
levels (the first test of Cook's proposal) [4], rotation of 
photon polarization [5], Landau-Zener tunneling [6], nu- 
clear spin isomers pp, level dynamics of individual ions 
[8], optical pumping [9], preservation of spin polariza- 
tion in gases [10 , quantum computing qubits undergo- 
ing decoherence [11] , Bose-Einstein condensates [12j [13] , 
optical systems [14] , NMR [15] , control of decay in opti- 
cal waveguides [16] and cavity quantum electrodynamics 
(CQED) [17 . Other experiments have also been pro- 
posed, involving neutron spin in a waveguide [18 and 
superconducting qubits [T9] . 

Remarkable applications of the QZE have been re- 
alized or proposed, such as the control of decoher- 
ence QTJ [20] , state purification [21] , implementation of 
quantum gates [22] and entanglement protection [23] . 
QZE can also inhibit entanglement between subsystems, 
making a quantum evolution semi-classical [23] . Other 
proposed applications consist in radiation absorption 
reduction and dosage reduction in neutron tomogra- 



phy [25], control of polarization [26] , and other general 
strategies to control decoherence [27] . 

In all these experiments or experimental proposals, re- 
peated projective measurements block the evolution of 
the quantum system in a non-degenerate eigenstate of 
the measured observable, so that the system is frozen by 
QZE in its initial state. However, more general phenom- 
ena can take place, for example when the measurement 
does not confine the system in a single state, but rather 
in a multidimensional (quantum Zeno) subspace of its 
Hilbert space. This gives rise to a quantum Zeno dy- 
namics (QZD) [28] : the system evolves in the projected 
subspace under the action of its (projected) Hamiltonian. 

No experiment has been performed so far to test the 
QZD. This would be important in view of possible ap- 
plications, for example in decoherence and quantum con- 
trol. We proposed in [29 a possible implementation of 
QZD in a CQED experiment. In this proposal, the field 
in the cavity undergoes a QZD under the joint action of 
a coherent source coupled to the mode (responsible for 
the Hamiltonian coherent evolution) and of a repeated 
photon-number selective measurement or unitary evolu- 
tion. This process is based on the spectroscopic interro- 
gation of the dressed levels of a single atom coupled to 
the cavity mode. These repeated operations create two 
orthogonal subspaces in the field's Hilbert space, with 
photon numbers larger or smaller than a chosen value s. 
QZD takes place in one of these subspaces. 

We also proposed in [29] that this procedure could lead 
to interesting methods towards the synthesis and manip- 
ulation of non-classical states. In this Article, we explore 
these ideas even further and detail some subtle mech- 
anisms involved in the dynamical evolution inside the 
Zeno subspace. We show, in particular, that the quan- 
tum Zeno dynamics can be used to produce mesoscopic 
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field state superpositions (MFSS), quantum superposi- 
tions of coherent components with different amplitudes. 
Such highly non-classical states are quite interesting for 
explorations of the quantum-classical boundary [30] . 

We start by introducing notations and by sketching 
the main ideas in Sec. [TTJ We explore the mechanisms of 
the confined dynamics and introduce the key idea of the 
'exclusion circle' in phase space in Sec. |III| The notion of 
phase space tweezers and scenarii of state manipulation 
are analyzed in Sec. |IV| Finally, we look at interesting 
perspectives on state synthesis in Sec. [V] We further dis- 
cuss practical implementation, that can be realized with 
a state-of-the-art apparatus in Sec. [VlJ where we also 
compare orders of magnitude and perform a few realis- 
tic simulations. Conclusion and perspectives are given in 
Sec. EEH 



II. GENERAL PRINCIPLES 

In this Section, we describe the principle of a quan- 
tum Zeno dynamics experiment in the cavity quantum 
electrodynamics context. The first Subsection ( II A[ ) ex- 
poses the general principle of the method and introduces 
useful notations. The method could be used in a vari- 
ety of experimental settings, and particularly in circuit 
QED [31J. However, for the sake of definiteness, we will 
discuss it in the framework of a microwave CQED experi- 
ment in construction at Ecole Normale superieure (ENS) 
involving circular Rydberg atoms and superconducting 
millimeter- wave cavities. We discuss the general features 
of this experiment in Subsection |IIB| We then describe 
how QZD may be implemented in this framework using 
repeated photon- number selective measu remen ts (IIC) 
or photon- number selective unitary kicks (II D). 



surement. The action of this measurement is represented 
by the projectors P^ [32], corresponding to the obtained 
result \i (J^ P M = 1). If the initial state is contained in 
the eigenspace associated to /io, the measurement gives 
almost certainly /i for each sequence, in the large TV and 
short r limit. The evolution is then confined in the Zeno 
subspace defined by P^ , which is in general multidimen- 
sional. The evolution in this subspace reads: 



U { P N \t) = [P, U{t/N)\ 



N 



— » e 



iHzt/hp 



for N — >• oo, where 



(1) 



(2) 



is the Zeno Hamiltonian. 

In the second procedure, the system undergoes a stro- 
boscope evolution, alternating short unitary evolution 
steps, governed by E/(r), with instantaneous unitary 
'kicks' Uk- The succession of N steps yields the unitary: 



Ug\t) = [U K U(t/N)] N ~ U%e- iH ^ h 



for N — >• oo, where 



Hz = J2 P » HP » 



(3) 



(4) 



Generalities and notation 



the P M s being the (multidimensional) eigenprojections of 

U K (U K P„ = e iX »Pjm- 

We observe that, by suitably choosing P Mo or Uk in 
Eqs. and ([3| respectively, one can modify the sys- 
tem evolution by tailoring the QZD, leading to possi- 
ble remarkable applications. We shall analyze here both 
schemes and discuss the experimental feasibility of the 
procedure Q-Q, related to the so-called 'bang-bang' 
control [33] used in NMR manipulation techniques [34] . 
The related mathematical framework is familiar in the 
context of quantum chaos [35 . 



A QZD can be achieved either by repeated (possibly 
unread) measurements of an observable with degener- 
ate eigenvalues, leading to a non-unitary evolution, or 
by repeated actions of a Hamiltonian kick with multi- 
dimensional eigenspaces, always leading to a global uni- 
tary evolution. The two procedures can be shown to be 
equivalent in the N — >• oo limit, where N is the num- 
ber of operations in a finite time interval t [20]. For N 
finite, differences can appear between the unitary and 
non-unitary procedures. Both measurements and kicks 
are supposed to take place 'instantaneously', namely on 
a timescale that is the shortest one in the problem at 
hand. We will discuss here both procedures before focus- 
ing on the latter, whose implementation in CQED turns 
out to be the easiest. 

The first procedure consists in N repetitions of a se- 
quence involving the evolution under the action of a 
Hamiltonian H for a time r = t/iV, generating the uni- 
tary U(r) = exp(—iHr/h), followed by a projective mea- 



B. A Cavity-QED setup 

Our proposal for QZD implementation [29] is based on 
the photon-number selective spectroscopic interrogation 
of the dressed levels for a single atom coupled to a high- 
quality cavity. In the ENS experiments, a very high-Q su- 
perconducting millimeter-wave cavity is strongly coupled 
to long-lived circular Rydberg states. The long lifetimes 
of both systems are ideal for the realization of experi- 
ments on fundamental quantum effects [30] . 

In all experiments realized so far, the atoms were cross- 
ing the centimeter-sized cavity mode at thermal veloci- 
ties (~ 250 m/s). The atom-cavity interaction time is 
thus in the few tens of /is range. It is long enough to 
result in an atom-cavity entanglement and short enough 
so that atoms crossing successively the cavity carry a 
large flux of information about the field state. This in- 
formation can be used for the implementation of ideal 
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FIG 



1. Scheme of the planned ENS cavity QED setup. 



quantum measurements [36] or for quantum feedback ex- 
periments [37]. However, this short interaction time is 
not compatible with a photon-number selective interro- 
gation of the dressed level structure at the heart of our 
QZD proposal. 

The ENS group is thus developing a new experiment 
with slow Rydberg atoms interacting for a long time 
with the cavity mode. Its scheme is represented on fig- 
ure [l] The Fabry- Perot cavity C [38 is made up of two 
superconducting mirrors facing each other (only one is 
shown in figure [l] for the sake of clarity). It sustains a 
non-degenerate Gaussian mode at a frequency close to 
51.1 GHz (6 mm wavelength). The mode has a Gaus- 
sian standing-wave envelope, with a waist w = 6 mm. 
Field energy damping times T c up to 130 ms have been 
reached by cooling the mirrors down to 0.8 K. At this 
temperature, the residual blackbody field corresponds to 
n t h = 0.05 photons in the mode on the average. 

The cavity is resonant with the transition between the 
two circular Rydberg levels e and g, with principal quan- 
tum numbers 51 and 50 respectively. These levels have 
a lifetime of the order of 30 ms, much longer than the 
typical atom-cavity interaction times considered in this 
Article (up to a few ms). Atomic relaxation thus plays a 
negligible role. 

The atoms are prepared by laser and radio-frequency 
excitation [30 out of a slow vertical atomic beam crossing 
the cavity in an atomic fountain arrangement. A Raman 
velocity selection performed on the slow beam emanating 
from a 2D-MOT source placed under the cavity makes it 
possible to selectively address atoms that are near the 
turning point of their ballistic trajectory at the cavity 
center. They thus reside in the mode's waist for a time 
of the order of 10 ms, limited only by their free fall. 

Excitation lasers are focused in C, delimiting a small 
volume. The initial position of the atoms is thus well- 
known. The time required for the atomic preparation is 
short (about 50 /is). It is important to note that this 
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FIG. 2. Dressed states of the atom-cavity system. The arrow 
indicates the photon-number selective transition addressed by 
S' for s = l. 



preparation does not involve any field close to resonance 
with the cavity mode. It can thus be performed without 
affecting the field quantum state. 

At the end of their interaction with the field, the atoms 
can be detected by the field ionization method inside the 
cavity itself. They are ionized by a field applied across 
eight electrodes circling the cavity and the resulting ions 
are routed towards a detector, which produces a macro- 
scopic signal. The method is state-selective, since the 
ionizing field depends upon the principal quantum num- 
ber. A simpler scheme can be used to perform an unread 
detection of the atoms, by merely ionizing them with a 
field applied directly across the cavity mirrors. Note that 
the centimeter-sized gaps between the ionizing electrodes 
enable us to couple millimeter-wave sources to the atoms 
or to the cavity mode (through its residual diffraction 
loss channels). 

In the following Subsections, we show how this ba- 
sic setup can be used for implementing the two QZD 
modes introduced in Section II A[ involving either re- 
peated photon-number selective measurements or re- 
peated photon-number selective fast unitary evolutions. 



C. QZD by repeated measurements 

The coherent evolution of the field in C is produced 
by a classical source S resonantly coupled with C [30] 
(Fig. [T]) . This evolution is described by the Hamiltonian 
(we use an interaction representation eliminating the field 
phase rotation at cavity frequency): 



H = aa^ 



a a 



(5) 



where a is the source amplitude and a (a^) the photon 
annihilation (creation) operator. If this evolution pro- 
ceeds undisturbed for a time interval t, a coherent state 
with an amplitude £ will 'accumulate' in the cavity, under 
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the action of the unitary displacement operator: 

U(t) = L>(0 = exp(f a 1 " - £*a) , f = -icrt/fi. (6) 

We now periodically interrupt this evolution over a to- 
tal time t by N measurements performed at very short 
time intervals r = t/N such that = | — iar/H\ <C 1. 
Each measurement involves a new atom prepared initially 
in the circular level h, with principal quantum number 49. 
Microwave pulses produced by the source S' probe the 
transition from h to g, at a frequency close to 54.3 GHz. 
Note that this transition is widely out of resonance from 
C. It can thus be probed without altering the mode state. 
Moreover, level h is impervious to the cavity field. 

Level g instead is strongly coupled to the cavity mode. 
In the resonant case, the atom-cavity Hamiltonian in the 
interaction picture reads: 



V 



(\e)(g\a + \g)(e\a 



(7) 



where is the vacuum Rabi frequency (Q/2ir = 50 kHz). 
The atom-cavity Hamiltonian eigenstates are the dressed 
states: 



M), |±,7i) 



1 
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\e,n-l)±\g,n)) , n > 1 , (8) 



where the former (latter) entry in each ket refers to the 
atom (cavity mode). The splitting between the dressed 
states |±,n) is hQy/n. 

The pulse sent by the source S' thus actually probes 
the transition between the level |ft,n) (whose energy is 
independent of the atom-cavity coupling) and the dressed 
states |zb, n). The level structure is shown in Fig. 2j The 
frequency of the \h,n) — >• transition depends upon 

the photon number n. 

Let us chose a specific photon number s > 1. The 
source S' is tuned to perform a (j) = tt Rabi pulse on the 
\h,s) — >• |+, s) transition. It is detuned from the bare 
h —> g transition frequency by In principle, we 

can chose the amplitude and duration St of this interro- 
gation pulse so that it has no appreciable effect on the 
transition between \h,s) and | — , s) or between \h,n) and 
|+, n) with n ^ s. This requires 1/St <C fl\ V s ± 1 — y^l- 
A long enough atom-cavity interaction time is thus essen- 
tial for the selective addressing of a single dressed atom 
transition. 

Finally, the source S' ideally performs the transforma- 
tions: 



and 



U 8 \h,s) = -i\+,s), 
U 3 \+,s) = -i\h,s) 

Us\-s) = \-s). 



(9) 



(10) 



If the cavity contains a number of photons n different 
from 5, then: 

U 8 \h,n) = \h,n), 

U a \±,n) = \±,n), (n^s) . (11) 



In conclusion. 

U a = 
P± = 



-i(\h,s)(+,s\ + \+,s)(h,s\) + P±, 
l-\h,s)(h,s\-\+,s)(+,s\ . (12) 



This is a unitary process: U s Ul = U^U S = t . 

We now examine the global evolution. Assume that 
the cavity is initially in its ground state and the atom in 
h. The joint atom-cavity state is \h,0). After the first 
time interval r, it becomes: 



e -«^/»|fc,0) = X>(r)|/i,n> , 



(13) 



where H is given by Eq. Note that H does not 

involve any atomic operator. Therefore, h is not affected 
by the coherent cavity evolution. The coefficients c n (r) 
are those of a coherent state with a small amplitude (3 = 
—iar/h. After this 'free' evolution for a short time r, the 
atom undergoes the tt Rabi pulse driven by S' (TT2J) : 



U s e~ iHT ' h \h,Q) = J2c n (T)\h,n)-ic s (T)\+,s) . (14) 

At this point, the atom is detected inside the cavity and 
its state recorded. Since \f3\ is very small, the probability 
for having s photons or more is small. With a large 
probability, the atom is thus found in h. This is our 
measurement: it makes sure that the number of photons 
in the cavity is not s. The cavity field is accordingly 
almost always projected onto: 



{h\U s e- iHT ' h \h) = Y J C n (T)\n) , 



(15) 



within a trivial normalization factor. One can thus sum- 
marize the action of U s followed by the measurement of 
\h) by the projection: 



P=t-\s){s\ 
acting on the field state alone, since 

Pe- iH ^\0} = J2cn(r)\n) 



(16) 



(17) 



[in terms of operators, we get from (12) that, in the pho- 
ton Hilbert space, (h\U s \h) = t — \s)(s\]. 

The Zeno procedure consists in the alternating evolu- 
tion under the action of the free Hamiltonian ^ and the 
projection (16): 



U ( p N) (t)= (Pe- iHr ' h ) N , t = t/N 



(18) 



which has to be understood as an evolution of the cavity 
field only. When N is large, one gets Up*\t) — >> Uz(i), 
where 



U z = e - iHzt/h P 



(19) 
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is the QZD generated by the the Zeno Hamiltonian 
Note that 

P = P<s + P>s , (20) 

where P <S (P >S ) is the projection onto the photon num- 
ber states with less (more) than s photons. Since H can 
create or annihilate only one photon at a time, one has 
P <S HP >S = 0, whence 



H z = P <S HP <S + P >a HP : 



>s 



H. 



<s 



(21) 



Here H <s is the restriction of the Hamiltonian H to the 
photon subspace H <s = P <s / H 1 spanned by the photon 
number states |0), . . . , \s — 1), and H >s the restriction to 
the subspace % >s containing more than s photons. 

Under the QZD, field states restricted to % <s and % >s 
remain confined in these subspaces, \s) realizing a hard 
'wall' between them. Strictly speaking, PH = H <S ^H >S 
forms a single Zeno subspace, within which evolution is 
coherent. For example, the coherence of the state (|n) + 
|p))/V% with n < s and p > s is fully preserved under 
the Zeno dynamics. However, transitions between H <s 
and 1~L >S are forbidden, due to the form of the interaction 
Hamiltonian. If the initial state is contained in only one 
of the two sectors, % <s or H> s , it will be confined to it. 
In the following, we shall focus on this situation. Note 
that, if C is initially in the vacuum state, with 5 = 1, the 
system remains inside H<i, i.e. in |0), and we recover the 
QZE E|. 

Of course, for a finite r (hence, a finite j3) there is 
at each step a small but finite probability for finding s 
photons in the cavity and, hence, the atom in another 
state than h. In this case, the cavity field is projected 
onto the s photon Fock state, and the Zeno procedure is 
abruptly interrupted. The probability of occurrence of 
such an event goes to zero when r is close to zero. For 
practical purposes, r should be chosen small enough to 
make this event unlikely in the planned duration of the 
experiment. 

Note also that, as long as the QZD is not interrupted 
by such an event, the atom is always found in h. It is 
thus not mandatory to actually record the atomic state at 
the detection time. Merely tracing over the final atomic 
state leads to the same results for the QZD of the field (in 
practical terms this means that the atom can be simply 
ionized by a field applied across the cavity mirrors, before 
being replaced by a new one for the next QZD step). 
In fact, as shown in the next Subsection, it is not even 
necessary to detect the atom at all. 



appreciably affect |ft,n) with n ^ s. This yields the 
transformation: 



U a \h,n) = (-l) s "-\h,n) , 



(22) 



and U s = 1 on all the other states. The atom always 
ends up in h, while the field experiences the selective 
kick Uk — U s with: 



U s = t-2\s)(s\ , 

UsUt = UtUs 



(23) 



1 . 



Such a photon-number dependent Rabi pulse [39 was 
used with s = 1 for a single-photon QND detection [40] 
and for a CNOT gate in CQED [41]. The evolution @ 
reads: 



U { P(t) = [U s U(t)} n - t/fe- iHz */ ft 
where t = Nr. The Zeno Hamiltonian Q is: 



M 7 



P_ = \s)(s 



>S j 



(24) 



(25) 
(26) 



which also satisfies Eq. (|21|. Once again, there is a hard 
wall at n = s, preventing transitions between % <s and 



E. Interrogation by a generic Rabi pulse 

We have seen that a QZD can be obtain ed both when 
S' drives a tt Rabi pulse, as in Sec. pICl or a 2n Rabi 
pulse, as in Sec. |IID| We show here that a generic pulse 
with an arbitrary Rabi angle (j) yields essentially the same 
physical effects. 

For a generic Rabi pulse, S' performs a unitary kick 
acting on the atom-cavity system, which mixes \h,s) 
with | +, s) and would create atom-field entanglement if 
C would contain s photons. The corresponding unitary 
operator reads: 



U s = exp 



-i^(\h,s)(+,s\ + \+,s)(h,s\) 



(27) 



For cj) = 7r, it reduces to the unitary of Sec. II C[ while it 
reduces to that of Sec. |IID| for <p = 2tt. The diagonaliza- 
tion of U s leads to: 



Us 



(28) 



D. QZD by repeated unitary kicks 

The same Zeno dynamics can be implemented by mak- 
ing use of a single atom, without any detection. Now, S f 
is tuned to perform a 2tt Rabi pulse on the s) —> |+, s) 
transition. As before, the pulse amplitude is weak enough 
(and its duration correspondingly long enough) not to 



with 

P± = \u ± )(u ± \, M= IM ^ + ' S) , (29) 

and Pj_ = 1 - P + - P-. In the large iV limit, for <j> ^ 0, 
the Zeno dynamics is generated by: 



H z = P+HP+ + P-HP- + P±HP ± . 



(30) 
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Since (u±\a\u±) = 0, the Zeno Hamiltonian reduces to 
H z = P ± HP ± . (31) 



The unitary (27) admits an invariant subspace of the 
range of the eigenprojection P± belonging to the eigen- 
value +1. Its projection is \h)(h\ ® (P <s + P >s ), the same 
as for a 2tt pulse. Starting from an atom in \h) and a field 
in 1~L <S or H> s , we obtain a QZD leaving the atom in \h) 
and the field in its initial subspace. 

Under perfect QZD, the cavity never contains s pho- 
tons and the atom and field are never entangled by the 
interrogation pulse. This discussion holds in principle for 
all non-zero values of (j) (0 < 4> < 2ir) in the N — »> oo limit. 
For finite values of TV, QZD is not properly achieved if <fi 
is very small, each kick operation being too close to 1. 
Numerical simulations, to be presented in Section |IIID[ 
fully confirm this qualitative argument. 



III. CONFINED DYNAMICS IN QZD 

We have shown that the QZD establishes a hard wall 
in the Hilbert space, corresponding to the Fock state \s). 
In qualitative terms, this hard wall can be viewed in the 
phase space (Fresnel plane) as an 'exclusion circle' (EC) 



with a radius y/s. In this Section, we examine the QZD 
starting with an initial coherent field located either in- 
side or outside the EC. This deceptively simple situation 
leads to the generation of a non-classical MFSS, quantum 
superposition of distinguishable mesoscopic states. 



A. Phase space picture 

Summarizing the main results of the preceding section, 
the Zeno dynamics consists in replacing the 'free' Hamil- 
tonian ([5| with the Zeno Hamiltonian Q: 



H z = J2 P » HP » = p >sHP >s + P <S HP < 



(32) 



For definiteness, let us assume that s = 4, identifying the 
two subspaces Hz = P< s T~i = span{|0), |1), |2), |3)} and 
H' z = span{|5), |6), |7), . . .}. In the photon- number states 
basis {|n)}, the only non- vanishing matrix elements of 
the annihilation and creation operators are: 



(n — l\a\n) = (n\aJ\n — 1) = \fn , 



(33) 



with n > 1. The matrix representations of H and Hz 
thus read: 



H = 





a 
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= Hy, 



(34) 



respectively. The two subspaces Hz and do not com- 
municate anymore in an ideal QZD situation. If the sys- 
tem is initially in Hz, it cannot make a transition to 1~L' Z 
and vice versa. Of course, this simple picture holds only 
in the limit of a true QZD, with an infinitesimally small 
time r between kicks and a very small displacement per 
step (3. We numerically examine the validity of this ap- 
proximation in Section [ill D| 

We simulated in [29] the QZD in cavity QED by apply- 
ing Eq. p4] ) [as already explained, Eq. ( p~8| ) would have 
been equivalent] with the unitary (|6| and by making use 
of Wigner's representation in phase space [42] . 

We summarize here the main results of our simulations. 
Figure [3] presents three sequences of 10 snapshots of the 
Wigner function W(^) separated by intervals of 5 steps, 
for s = 6 and /3 = 0.1. The Wigner function is defined 



by: 

w (0 = ^nj^- v\e\ Re ^ + yy 2ylm(/h dy , (35) 

g being the density matrix of the photon field, obtained 
by tracing out the atomic variable. 

In Fig. ^ a) the field is initially in its vacuum state 
|0) G %<6- Its amplitude increases along the real axis 
(free dynamics). When this amplitude reaches ~ 2, be- 
tween 15 and 20 steps, the coherent state 'collides' with 
the EC of radius a/6 (dashed line in Fig. [3]). The field am- 
plitude stops growing and undergoes a very fast 7r phase 
shift between steps 20 and 30. At step 25, the field is 
in a MFSS, quantum superposition of two components 
with opposite phases. The fringes inside the EC are the 
signature of the quantum coherence. At step 35, the field 
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FIG. 3. (a) QZD dynamics in 1~L<q. Ten snapshots of the field Wigner function W(£) obtained after a number of steps indicated 
above each frame. The cavity is initially in its vacuum state, s — 6 and f3 = 0.1. The EC is plotted as a dashed line, (b) QZD 
dynamics in / H > q. Same as (a) with an initial a = — 5 amplitude, (c) Same as (b), with an initial amplitude a = — 4 + In 
(b) and (c) the successive frames correspond to the same step numbers as in (a). From [29] . 



state is nearly coherent again with an amplitude close to 
—2. It then resumes its motion from left to right along 
the real axis, going through |0) around step 45 and head- 
ing towards its next collision with the EC. The long-term 



dynamics will be discussed in Section IIIC 



QZD in H>6 is illustrated in Fig. [3^b), with snapshots 
of the field Wigner function for s = 6 and an initial coher- 
ent state \a = —5). The field collides with the EC after 
20 steps. It undergoes a QZD-induced 7r phase shift be- 
ing, after 25 steps, in a MFSS. After 30 steps, the state 
is again nearly coherent with a positive amplitude and 
resumes its motion along the real axis. After 45 steps, 
its amplitude is slightly larger than 4.5. It would be —0.5 
in the case of free dynamics. 

Finally, in Fig.[3^c), the field state collides tangentially 
with the EC. The parts of the Wigner function that come 
closest to the EC propagate faster than the others. The 
state is distorted and squeezed (albeit by a moderate 
amount) along one direction. 



Phase inversion mechanism 



We now show that the main feature of the evolution of 
the Wigner function, the fast phase shift during the colli- 
sion with the EC, can be understood via a semi-classical 
argument. Let us set a = i/y/2 in ^ for simplicity. We 
get that the Hamiltonian, 



is simply the momentum operator. Thus, by the spectral 
theorem, 



s-1 



H <s =P <sP P <s = \n)(n\p\n')(n'\ 

n,n'=0 

= X[0,J5 a _i](#h.o.) P X[0,£ s _i](#h.o.) , (37) 



where 



ffh.o. = ^(X 2 + P 2 ) 



(38) 



is the harmonic oscillator quantum Hamiltonian (with 
m = uo = 1), whose energy is forced to be less than 
E s -i = h(s — 1/2) by the characteristic function x 
[Xa(x) = 1 if x G A,0 otherwise]. 

For large quantum numbers, we can approximate H <s 
by its classical limit, which reads 



h(x,p) =PX[o,R](r) 



(39) 



where r = ^/ x 2 + p 2 and R = ^/2E s -i. This Hamilto- 
nian describes the motion of an ultra-relativistic particle 
(energy proportional to momentum) confined in phase 
space by a hard wall at r = R (non-holonomic con- 
straint). 

The Hamilton equations of motion are: 



x = X[o,R](r) -p6 R (r)- , 
x 

p = p S R (r)- , 



(40) 
(41) 



where Sr is the Dirac delta function at R. If the particle 
is not on the EC, this yields the solution 



H = i(a-a t )/V2 =p 



(36) 



p = Po , X = X + t , 



(42) 
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FIG. 4. Vector field in classical phase space, (a) Motion inside 
the EC; (b) Motion outside the EC. 



xq and po being the initial position and momentum, re- 
spectively. The particle is thus proceeding at a constant 
velocity along the x axis. When it hits the EC, the evo- 
lution is dominated by the singular contributions in (40)- 
(41) through the vector field: 



X(x,p) = - S R (r) 



(43) 



that yields a motion along the circle at a constant speed 
(infinitely large in the limit of an infinitely sharp confine- 
ment inside the EC). The particle reappears on the other 
side of the EC (with the same momentum) and resumes 
its motion along the x axis at a constant velocity. These 
trajectories are qualitatively sketched in Fig. [4ja). 

A cloud of such particles would thus evolve essentially 
as the field Wigner function inside the EC. This explains 
the 'phase inversion mechanism' of Fig.[3^a): the Wigner 
function hits the right hand side of the EC and almost in- 
stantaneously reappears on the left hand side. Of course, 
the transient creation of an MFSS involving a quantum 
superposition of two large fields with opposite phases and 
the appearance of an interference pattern inside the EC 
(Figure |3j frame 25) cannot be accounted for in this clas- 
sical picture. 

When the particle is initially outside the EC, the evo- 
lution is generated by the Hamiltonian: 



h(x,p) =p(l-X{o,R](r)) 



(44) 



and the conclusions are identical. When the particle hits 
the EC, it moves very quickly to the other side [Fig.Qb)]. 
This explains the fast motion of the components of the 
Wigner function that come closer to the EC, in Figs.[3jb) 
and (c). 



C. Long-term evolution 

We now analyze the long-term evolution of the field en- 
ergy when the state is initially inside the EC. In this case, 
only a finite set of Bohr frequencies is involved in the evo- 
lution, which is thus expected to be quasi-periodic [30] . 



State distortions, however, eventually accumulate and 
damp the oscillations of the field amplitude. This phe- 
nomenon was numerically investigated in [29 and will 
now be analyzed in greater detail. 

Without loss of generality, one can consider a real and 
write: 

H = a(o) +a) = a^y] Vn + 1 (\n)(n + 1| + \n + \){n\) . 



n>0 



(45) 

Indeed, Hamiltonians ^ and (45) are unitarily equiva- 
lent via U ((f) = e^ ata , ip being the phase of a in Eq. 
They have thus the same spectrum and generate the same 
dynamics. 

Let us look first at rather small values of s. For 8 = 4, 
all properties of the Zeno dynamics in Hz depend on 
those of the matrix 



H <s = P <S HP <S = a 



<s 



<s 



/0 1 \ 

1 y/2 

y/2 V3 

\0 a/3 / 



(46) 



that has four non-degenerate eigenvalues ±aA + = 

and ±aA_ = d=a \/3- a/6, so that the Zeno 
dynamics of a generic observable is a quasi-periodic mo- 
tion on a four-dimensional torus. We get: 

|(0|e-^<^|0)| 2 



Po(t) 

P2(t) 

Ps(t) 



((A+ - 1) cosujX-t + (1 - X 2 _) cosuj\+t) 



1 

24 

1 2 

— (A_|_ sinco>A_£ + A_ sinco>A + £) 

1 2 

— (cos uj\-t - coscjA+t) , 

1 2 

— (A+ sinujX-t — A_ sina;A + £) 



(47) 
(48) 

(49) 

(50) 



where uo = a/H. Thus, the population of the number 
states are quasi-periodic. We have: 



1 



6 



(51) 



so that A + /A_ ~ 22/7 up to a few per mil. Therefore, 
the populations almost return to their initial value after 
a time 



22- 



2tt 
ujX-l. 



2tt 
uj\- 



7- 



59.2 



(52) 



For larger values of 5, the calculations become analyt- 
ically unmanageable. However, the main conclusions re- 
main valid and the features of the evolution qualitatively 
identical. 

The average photon numbers as a function of uot for 
5 = 4 and s = 6 are displayed in Fig. [5| We observe 
at long times a quantum revival [30] at wt c± 59.2 and 
ut ^ 150, respectively. 
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FIG. 5. Number of photons as a function of uot = at/h. Upper 
panel: s = 4, lower panel: s — 6. Note in both cases the 
recurrence of the photon number oscillation at uot c± 59.2 and 
uot c± 150, respectively. 



D. Limits of QZD and applications 

We have so far assumed a perfect confinement inside 
the EC. It can be obtained only in the limit of vanishingly 
small displacements j3 at each step and for non- vanishing 
interrogation pulse Rabi angles <j). In a real experiment, 
the preparation and interrogation take a finite time and 
thus the displacement per step cannot be made arbitrar- 
ily small. We have explored the corresponding limits of 
the QZD by extensive numerical simulations. We focus 
here on the typical example of a dynamics inside the EC 
with 5 = 6 [Figure |3| a)]. 

The calculations have been performed using the quan- 
tum optics package for MATLAB [42 . The field Hilbert 
space is truncated to the first 60 Fock states. The initial 
cavity state is the vacuum. Each step involves a trans- 
lation by an amplitude f3 (chosen real positive without 
loss of generality). It is followed by a Rabi interrogation 
pulse with an angle (j) on the 6) — >• |+,6) transition 
for an atom remaining motionless at cavity center. No 
atomic detection is performed and the same atom is used 
for all elementary steps. No other experimental imper- 
fections (cavity relaxation, finite selectivity of the Rabi 
pulse, etc.) are taken into account. They will be dis- 
cussed in Section IVTl 

To assess the quality of the confinement, we compute 
the evolution for a number of steps N = I[2y/6/ /3] (where 
/ stands for the integer part), corresponding to the first 
return of the field state close to the vacuum. We compute 




FIG. 6. Fidelity of confinement for s — 6 versus the transla- 
tion per step f3 and the interrogation pulse Rabi angle <j) in 
radians. The contour lines mark the 95% and 98% levels 

the fidelity 

F = Tr( 6Qp ) (53) 

of the final field state g with respect to the reference state 
g p obtained after the same number of steps in an ideal 
QZD (evolution in an Hilbert space strictly limited to the 
first 6 Fock states). 

Figure [6] presents F as a function of f3 from 0.05 to 
1 (0 is obviously excluded) and of <j) from to 2tt. Not 
surprisingly, F is close to one for small j3 and large (j). 
It is nearly zero when (j) approaches zero or (3 one. The 
contours correspond to 95% and 98% fidelities respec- 
tively. As shown in Section |VI[ large <ft values can be 
easily implemented in a short time interval. For a 2tt 
interrogation pulse, we can achieve an excellent confine- 
ment fidelity (98%) with a translation per step as large 
as /3 = 0.4. The QZD is thus a quite robust mechanism. 
This is promising for practical applications. 

When F is close to 50%, an interesting situation arises. 
We observe numerically that, when the collision of the 
moving coherent state with the EC occurs (after a num- 
ber of steps around AT/2), part of the Wigner function 
is transmitted through the barrier. Another part under- 
goes the phase inversion mechanism and is rejected to 
the other side of the EC. In the following steps, these 
two components evolve separately in the subspaces H<6 
and H>6- The outer one moves further along the posi- 
tive real axis and the inner one returns close to the origin. 
Finally we are left with two nearly coherent components 
centered at the origin and at an amplitude 2\/6. 

The final trace over the atomic state does not erase the 
coherence between these two components, showing that 
the atom is not strongly entangled with the field even 
though the \s) Fock state has been transiently populated 
in the process. We are thus left with a quite mesoscopic 
field state superposition. The evolution of the Wigner 
function corresponding to f3 = 0.345 and <j) = 3.03 is pre- 
sented on figure [7](a). The coherence between the two 
components in the final MFSS is manifest with the pres- 
ence of the characteristic interference fringes [figure [7^b)], 
even though the contrast of these fringes is not maximal. 
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FIG. 7. Generation of a cat by a semi-transparent EC 
(f3 = 0.345 and <f> — 3.03). (a) 3 snapshots of the field Wigner 
function W(£). The corresponding number of steps are indi- 
cated above the frames, (b) Final Wigner function after 14 
steps. The negative parts are conspicuous. 




FIG. 8. Transmission of the EC as a function of the QZD 
parameters. The solid green line follows the 50% transmission 
level. The conditions used for figure[7]are marked by the white 
dot. 

component on a partially-transparent EC. The relative 
weights of these components can be adjusted by fine- 
tuning, during each collision, the incremental step f3 and 
the interrogation pulse angle <ft. By changing the phase 
of j3 from one collision to the next, a wide variety of 
multi-components superpositions can be produced. 



We have systematically studied the generation of such 
MFSS by an imperfect QZD. For each value of /3 and 0, 
we compare the final cavity state to a superposition of 
two coherent states : 

|MFSS) = w <s \a <s ) + w >s e ie \a >s ) , (54) 

with real amplitudes a <s and a >s for the confined and 
transmitted parts respectively (a <s being close to zero 
and a >s close to 2y/6). These amplitudes, the two real 
mixture coefficients w <s and w >Sl and the relative quan- 
tum phase are fitted to optimize the fidelity of this ref- 
erence state with respect to the final state in the cavity. 
We find that the relative phase is always very close to 
7r. For the conditions of figure [7[ the fidelity is 75%. It is 
limited in particular by a residual spurious entanglement 
between the atom and the field. 

We define the EC 'transparency' T as: 

T = w 2 >s /(w 2 <s + w 2 >s ) , (55) 

weight of the transmitted component in the final MFSS. 
Figure [8] presents T versus the QZD parameters /? and 0. 
The 50% level is indicated by the green line and the con- 
ditions of figure [7] by the white dot (transmission 44%). 
We observe that MFSS can be generated in a large range 
of operating parameters. Note that, for very large j3 val- 
ues, the state transmitted through the barrier can be 
notably distorted. The fidelity with respect to an ideal 
MFSS is then rather low. 

This MFSS generation method can be straightfor- 
wardly generalized to superpositions of more than two co- 
herent components by repeated collisions of the trapped 



E. QZD in a translated EC 

The QZD proceeds in an EC centered at the origin 
of phase space. It can be straightforwardly generalized 
to an EC centered at an arbitrary point in phase space. 
Before the interrogation pulse, we perform with the help 
of the source S a displacement of the field by a (possibly 
large) amplitude —7. After the interrogation pulse, we 
translate back the whole phase space by the amplitude 7. 
Qualitatively, we block the evolution in an EC centered 
at the origin for a field state globally translated by the 
amplitude —7. This is clearly equivalent to blocking the 
evolution in an EC centered at the point 7 in phase space. 

In more precise terms, the kick operator Uk is changed 
by the two translations from U s into 

U,{i) = D( 7 )U s D(- 7 ) . (56) 

After p steps, the global evolution operator is 

U z (s, 1 ,p) = [U s ( 1 )D(f3)}r , (57) 

which can be expressed, using displacement operator 
commutation relations, as: 

U z (s,j,p) =D( 7 )C/z(s,0,p)D(-7)exp[2ipIm(/? 7 *)] . 

(58) 

Up to a topological phase, the state after p steps is equiv- 
alently obtained by first displacing the field by —7, then 
performing p QZD steps in an EC centered at origin and 
finally displacing back the field by 7. 

The s = 1 case is particularly interesting in this con- 
text. The QZD blocks the unique coherent state I7) at 
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a fixed point in phase space, while all other parts of the 
phase space can be moved by the action of displacement 
operators. This ability to operate separately on differ- 
ent regions of the phase space will be instrumental in the 
next Section. 



IV. PHASE SPACE TWEEZERS 

We proposed in [29] to use an s — 1 EC as phase-space 
tweezers. Let us assume that the initial state of the cavity 
field is made up of a superposition of non-overlapping 
coherent components, prepared for instance by using in 
a first stage of the experiment a semi-transparent EC. We 
can use an s = 1 exclusion circle to block one of these 
components, with an initial amplitude 70. We assume 
here that, besides the displacements used to generate the 
off-center EC, there is no other source of evolution of 
the field. Now, we change at each step of this new QZD 
dynamics the center of the exclusion circle, from 70 to 
71, . . ., to 7AT. 

Provided the difference between two successive posi- 
tions of the EC, |7i+i — 7i|, is always much smaller than 
one, the coherent component trapped in the EC will adi- 
abatically follow the motion of its center. The coher- 
ent state amplitude will thus be changed from 70 to 7 at, 
while all other components of the initial state remain un- 
changed. 

The movable EC operates in phase space as the opti- 
cal tweezers which are now routinely used to move mi- 
croscopic objects. In analogy, we coined the term 'phase 
space tweezers' for this operation [29] . 

Obviously, ideal operation of the phase space tweez- 
ers requires an infinitely small increment of the EC po- 
sition at each step, hardly compatible with a practical 
implementation. We have thus studied the quality of the 
tweezers operation with respect to the interrogation pulse 
characteristics and to the 'velocity' of the EC motion. 

We compute the final field state for a tweezers action, 
taking initially the vacuum state and pulling it away. 
As in Section |III D all sources of experimental imper- 
fections are neglected in this calculation. The exclusion 
circle center moves in N steps from zero to a real am- 
plitude 2^/6 (i.e. 24 photon field, this amplitude being 
chosen rather arbitrarily to coincide with the total field 
displacement used in Section HID). We compute the fi- 



nal state fidelity with respect to a coherent state with 
an optimized amplitude. We observe in fact that, for 
small N values i.e. large EC displacements per step, the 
state wiggles slightly inside the EC during translation. 
The final amplitude might thus not be exactly equal to 
2^. The effect is quite small, the maximum amplitude 
difference being less that 0.5 for all data presented here 
Figure [9] presents the calculated final fidelity as a func- 
tion of the number of steps N (from 10 to 31) and of the 
interrogation pulse Rabi angle <p (from zero to 2ir). The 
95% and 98% levels are indicated by the blue and red 
contour lines. For the 2tt interrogation pulse, the fidelity 
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FIG. 9. Tweezers operation fidelity as a function of N and 
<j). The fidelity is computed with respect to an optimized 
coherent state. The two solid lines follow the 95% and 98% 
levels. 




FIG. 10. (a)-(b) Initial and final Wigner functions W(£) for 
a phase space tweezers operation. The first steps ECs are 
depicted as solid lines in (a) and dotted lines in (b), the final 
ECs as solid lines in (b). The arrows in (b) indicate the two 
EC centers trajectories. From |29| . 



is extremely large, more than 99% for N = 10, i.e. a 
motion of nearly 0.5 per step. Once again, the QZD op- 
eration is very robust and tweezers can be used to move 
quite rapidly a coherent component through the phase 
plane. 

Obviously the fidelity decreases quite rapidly with the 
interrogation angle <j). For <\> = 7r, N = 60 steps are 
required to achieve a 99% fidelity. When the EC step 
is too large, or the interrogation angle to low, the EC is 
slightly transparent at each step and leaks a little bit of 
the trapped state. The final state is stretched along the 
path of the EC, with no remarkable features. 

Phase space tweezers can be used to increase at will 
the distance between two components of an MFSS as 
The initial cavity state is \a) + 



illustrated in figure 10 
I — a) with a = 2. It is turned in 100 steps (50 for the 
motion of each component) into + | — a') with a' = bi. 
The final fidelity is 98.8% with respect to the expected 
cat. A wide variety of operations on non-classical fields 
can be envisioned with this concept. 

The tweezers operation allows us to tailor at will a 
pre-existing superposition of coherent states. It can be 
slightly modified to allow for the generation of this super- 
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position from the vacuum, as shown in the next Section. 

V. STATE SYNTHESIS 

We present in this Section a QZD-based method for 
the generation of a nearly arbitrary superposition of co- 
herent components, starting from the vacuum state. We 
proceed with an EC motion driving the vacuum to the 
first required coherent component. However, we perform 
a quantum superposition of the tweezers operation and 
no operation at all by casting the interrogation atom in 
a superposition of ft with an inactive state i (for instance 
the circular state with principal quantum number 52, ly- 
ing above e), which does not take part in the QZD pro- 
cess. The atom gets in the process entangled with a field 
involving a superposition of the vacuum with a moving 
coherent component. This process is then repeated for 
all the coherent components in the final state. 

Let us write the target state |\E^) as: 

m 

l*t>=£c J -| 7i > , (59) 
i=i 

superposition of m coherent states. We assume that all 
the 7jS with j ^ m have a negligible mutual overlap as 
well as with the vacuum state. Up to an irrelevant global 
phase, we can also assume c m real. 

We first create the component I71) out of the initial 
cavity vacuum state |0). The atom, initially in ft, is at 
rest at cavity center. We send on the atom, with the 
help of a microwave source 62, a narrow-band (soft) mi- 
crowave pulse resonant with the |ft, 0) — » \g, 0) transition, 
tuned to produce the state superposition a\ \g, 0) +61 1 ft,, 0) 
(the a\ and b\ coefficients will be determined later). 

A tweezers operation performed with the atom initially 
in g and an empty cavity leads to a partially transparent 
EC (the initial state has a component on | — ,0), which 
is not addressed by the interrogation pulses) and to a 
spreading in phase space. We must avoid this effect. Be- 
fore performing the tweezers action, we thus shelve level g 
in the fourth level i. We use for this purpose a millimeter- 
wave source £3 tuned to resonance with the two-photon 
transition between g and i at 2x49.6 GHz. The strong 
coupling of Rydberg atoms to millimeter-wave sources 
makes it possible to achieve a tt pulse on this transi- 
tion in a short time interval. Such a short (hard) pulse 
does not resolve the dressed level structures and performs 
the transition whatever the photon number in the cav- 
ity. Finally, we reach the quantum state superposition 
ai|i,0)+fti|ft,0). 

We then perform the tweezers action itself, using the 
interrogation source S' tuned for s — 1 and the trans- 
lation source S. The tweezers is active only if the atom 
is initially in state ft. The EC center evolving from 
to 71, we are finally left with the entangled atom-cavity 
state ai|z,0) + &i|ft,7i). We do not take into account 
here any topological phase that could affect the |ft, 71) 



part of the state if the trajectory through phase space 
was not a straight line. This phase could easily be taken 
into account with minor modifications of the algebraic 
expressions. A final hard — tt pulse on the i — » g transi- 
tion driven by 63 leads us to ai\g : 0) + 61 1 ft, 71). 

Since 71 is notably different from zero, a soft pulse on 
the |ft,0) — >• |g,0) transition driven by 62 addresses only 
the part of the atom-cavity state involving the vacuum. 
We tune this pulse to produce the state superposition 
ai{a 2 \g, 0) + 62 1 ft? 0)) + 61 1 ft, 71). We then shelve g in i 
with a hard pulse driven by 63 and perform a tweezers 
operation leading from the vacuum to the amplitude 72. 
We should take care that the EC never comes close to the 
71 component, which should be left unchanged. There is 
of course ample space in the phase plane to plan a conve- 
nient trajectory. Finally, we unshelve level z, leading to 
the state aia 2 \g,0) + aib 2 \h^ 2 ) + &i|ft>7i)> involving a 
superposition of three coherent components, two of them 
(71 and 72) being disentangled from the atomic state. 

Since again 72 is notably different from zero, we can 
selectively address with S 2 the g part of the state to split 
it in a coherent superposition. Iterating the process m 
times, we prepare finally the state: 

a t a 2 • • -a m _i|#,0) + aia 2 a 3 • • • a m _ 2 6 m -i |ft, 7m-i) + 
... + aift 2 |ft,72>+&i|^7i) • (60) 

A final tt pulse produced by S 2 on \g, 0) — » |ft, 0) casts the 
atom in ft with certainty and a final tweezers operation 
from to 7 m leaves the atom in |ft) and the cavity in the 
state: 

ai<2 2 • • -a m -l|7m) + «1«2«3 • • ■ CLm-2bm-lhm-l) 

... + aia 2 6 3 |7 3 ) +aifr 2 |72) + &1I71) • ( 61 ) 

We must now determine the intermediate coefficients 
di and bi so that the final state is \Sfr t ) [Eq. ([59])] . The 
si mplest ch oice is obtained by setting b\ = c\ and a\ = 
y/l — \bi\ 2 . This determines the value of b 2 and hence 
(within an irrelevant phase that we take to be zero), that 
of a 2 . We then get 63 and as and stepwise all the required 
coefficients. 

We have numerically simulated the procedure for the 
creation of a four-component MFSS: 

\ (|4> + |4i) + |3e^ 4 > + |0>) . (62) 

All tweezers actions are performed with a 0.1 amplitude 
increment and a (j) = 2ir interrogation pulse. The Wigner 
function of the resulting state is plotted in figure [TT] The 
fidelity with respect to the target state is 99%. 

This method opens many perspectives for the genera- 
tion of complex MFSS. The only restriction is that the 
final components should not overlap with each other and 
with the vacuum (if this is not the case, it is not possi- 
ble to manipulate one independently of the others with 
the tweezers). There is nevertheless a wide range of 
state superpositions that can be directly reached with 
this method. 
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FIG. 11. Final state Wigner function W(£) after the synthesis 
of a complex state superposition. See text for the conditions. 



VI. SIMULATIONS OF A REALISTIC 
EXPERIMENT 

We have up to now discussed the QZD in an ideal set- 
ting, assuming no atomic motion, no cavity relaxation 
and, more importantly, a perfect selectivity of the inter- 
rogation pulse. We now proceed to include a realistic 
description of these imperfections and to assess the qual- 
ity of the QZD in this context. 

First, atomic motion through the mode, at a low veloc- 
ity, is not a real problem. Provided the initial position of 
the atom (determined by the excitation lasers) and the 
atomic velocity are well known, the position of the atom 
is precisely known at any time during the sequence. The 
slow variation of the atom-cavity coupling can then be 
taken into account. We can for instance tune the inter- 
rogation pulse source to remain resonant on the selected 
dressed atom transition. 

We thus address in more details the two other issues. 
The main one is the interrogation pulse selectivity, dis- 
cussed in the next Subsection. The final Subsection is 
devoted to realistic simulations of a few key QZD exper- 
iments including cavity relaxation. 



A. Interrogation pulse optimization 

The interrogation pulse should address selectively the 
|ft, s) — >• |+, s) transition. For a motion inside the EC, 
or for a tweezers operation, this pulse should not af- 
fect any transition corresponding to a photon number 
n < s, and more particularly the transition between 
|ft, s — 1) and |+, s — 1), which is closest to resonance 
with the interrogation pulse. The frequency difference 
between the addressed and the spurious transitions is 
only S n = (Q/2)\y/s — y/n\. In particular, 5 s -i decreases 
with increasing s. 



For the sake of definiteness, we shall only consider two 
practical cases in the following, that of a tweezers oper- 
ation (s = 1) and that of a motion inside the EC with 
s = 6. In the latter case, the interrogation pulse should 
resolve a frequency splitting 5^/2tt = 5.3 kHz only. Rely- 
ing on a very long pulse duration to achieve this resolu- 
tion leads to unrealistically long experimental sequences 
in view of the finite cavity lifetime. A careful tailoring of 
the interrogation pulse is the only realistic solution. 



1. Square pulses 

The simplest procedure is to set S' to produce a square 
pulse with a duration t p , resonant with the addressed 
transition and performing a Rabi rotation by an angle 
n p 7r. In most cases, n p is set to two, but smaller values 
can be used (n p = 1 is appropriate to implement a semi- 
transparent EC - note that n p need not be an integer). 
To minimize unwanted transitions, we can chose the pulse 
duration so that the same pulse produces a p p 7r Rabi 
pulse on the non-resonant nearby transition (n = s — 1), 
where p p is an even integer (obviously larger than n p ). 
This condition sets a zero in the spectrum of the pulse 
at the precise frequency of the spurious transition. 

A simple algebra on Rabi rotations leads to a pulse 

duration t p = (n /5 s -i) ^Jp p — n p . For 5 = 1 and n p = 2, 

we get t p = 69 /as for p p = 4 and t p — 113 /is for p p = 6. 
In the more demanding s — 6 case, the pulse durations 
with the same settings are 324 and 530 /is respectively. 
All these durations are still much shorter than the atom- 
cavity interaction time scale. 

Numerical estimations of the influence of this square 
pulse on the complete dressed states structure confirm 
that the transfer rates on the non-resonant transitions 
are small, in the % range at most for the shortest pulses. 
However, we observe that the relatively strong pulse pro- 
duces an appreciable phase shift of the states \h,n) with 
n < s, due to the accumulated light shift effect. This 
is not a too severe problem for the tweezers operation, 
since this amounts finally to a global, predictable, phase 
shift on the displaced component. This phase shift could 
in principle be compensated for or taken into account in 
the state synthesis method. 

The influence of these light-shifts is, however, much 
more obnoxious for the QZD in an EC. Setting 5 = 6 
and p p = 4, we get a —0.86 rad shift for |ft, 5), —0.47 rad 
for | ft, 4). This phase shift is inversely proportional to S n . 
It is not a linear function of n and cannot be absorbed, 
as an index of refraction effect, in a mere redefinition 
of the cavity frequency. We have checked by numerical 
simulations that this phase shift destroys most of the 
QZD features. The EC remains an impenetrable barrier, 
but the state inside it is completely distorted, even far 
before it reaches the EC for the first time. 

Note that pulse shape optimization can reduce the spu- 
rious transfer rates but does not solve the light shifts 
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problems. A more sophisticated pulse sequence is manda- 
tory. 



2. Optimized composite pulses 

We propose thus to use a composite pulse sequence 
that leads to a nearly perfect cancellation of the light 
shifts. We discuss it in the important case of a 2tt inter- 
rogation pulse. The sequence is made up of three pulses: 

• a 7r square pulse on \h,s) — » carefully opti- 
mized as in the previous Subsection (n p = 1, p p = 2 
or 4). 

• a fast phase shift of the atomic levels alone, chang- 
ing \g) into —\g) and amounting to exchanging the 
|+, n) and | — , n) dressed states for all photon num- 
bers. 

• an optimized tt pulse on | — , s) — > \h,s), similar to 
the first one within an adjustable phase cp. 

In principle, the first pulse transfers all the population 
from \h,s) to |+, s). The central phase shift transforms 
|+, s) into | — , s). The final pulse transfers back the state 
into s), with a phase shift that can be adjusted by 
tuning the phase cp. For all other levels (n ^ s), there 
is almost no transfer out of \h, n) by the initial and final 
optimized pulses and the central operation has thus no 
effect. In this composite sequence, the two microwave 
pulses applied on the atom-cavity system have opposite 
detunings with respect to the spurious \h,n) —> |+,n) 
and | — , n) — » \h,n) transitions (n < s). One can thus ex- 
pect that the phase shifts due to the second pulse exactly 
compensate those produced by the first. 

The central phase shift operation could be performed 
by a hard non-resonant pulse coupling g to another level 
(i for instance). The accumulated light shift can be tuned 
for an exact tt phase shift of g, independent upon the pho- 
ton number in the cavity. Levels h and e, farther away 
from resonance with this dressing pulse, are not affected. 
A simpler solution is to use the differential Stark shift on 
the three levels e, g and h as in [43]. A short pulse of 
electric field applied across the cavity mirrors produces 
three different photon- number-independent phase shifts, 
tp e , tyg and cp h on these levels. Setting cp e — (p g — 7r, 
we are left within a global phase with the transforma- 
tions |e) — )► |e), \g) — >• —\g) and \h) — >> e z ®\h). For most 
QZD operations, the phase <I> acting on h is an irrelevant 
global phase factor. For the state synthesis, it is a well- 
known quantity that can be taken into account in the 
state preparation sequence. 

We have evaluated numerically the effect of this com- 
posite pulse. For s = l,n p = l,p p = 4, the total duration 
of the sequence is 155 /is. With ip = 2.75 rad, we perform 
the selective transformation 1) —> —\h,l). The resid- 
ual phase shift on 0) is only —3. 10 -5 rad. In the more 
demanding s = 6 case, we use n p = 1, p p = 2 for a total 



duration of 325 /is. The residual transfer rates for n < s 
are below 1.5% and the residual phase shifts lower than 
10 -4 rad. For the few tens of pulses in a typical QZD 
sequence, these imperfections have a negligible influence. 

The principle of this composite pulse can be extended 
to other <j> values and, in particular, to = 7r, useful 
for the realization of semi-transparent ECs. This pulse 
is made up of an n p = 1/2 pulse on \h,s) — » |+, s), fol- 
lowed by the atomic phase inversion and a n p — 1/2 
pulse on |— , s) — >• \h,s) with a phase cp. In the ideal 
case, this pulse combination results in the transforma- 
tions \h,s) —> | — , s); |— , s) —> |+,s) and |+, s) —> \h,s) 
(three applications of the transformation are necessary 
to return to \h,s)). Since the level \h,s) is nearly never 
populated in a successful QZD, this composite pulse is 
basically equivalent to a standard tt pulse. 

The composite pulse architecture achieves the required 
selectivity in a relatively short interrogation time. We 
use now these optimized pulses for the simulation of a 
few key QZD experiments. 



B. Simulation of key experiments 

The simulations include a realistic description of the 
composite interrogation pulse and of cavity relaxation. 
We use the best available cavity damping time T c = 
130 ms [38]. We should of course check that the total 
duration of the sequence remains in the ms time range, 
much shorter than the atomic free fall through the cavity 
mode. 

The periodic motion of a coherent component inside 
the EC is barely affected by cavity relaxation. Setting 
5 = 6, using a composite 2tt interrogation pulse with 
a translation per step f3 = 0.4, we get a fidelity after 
12 steps (corresponding to the return near the vacuum 
state) of 90% instead of 92% in the ideal case treated in 
Section |III| The total sequence duration is 3.9 ms. In 
fact, the field propagates most of the time inside the EC 
as a coherent state, nearly impervious to relaxation. It 
is only during the phase inversion, for a few steps, that 
a MFSS prone to decoherence is generated. 

We have also examined the creation of a MFSS by 
transmission trough a semi-transparent EC. With a com- 
posite 7r interrogation pulse, s = 6 and j3 = 0.33, we ob- 
tain in 5.4 ms a fidelity with respect to an ideal state of 
79%. This is promising to study the decoherence of this 
large MFSS (the square of the distance in phase space be- 
tween the two components, setting the decoherence time 
scale [30], is 24 photons). 

As a more striking example, we have simulated the 
generation of a three-component MFSS by two collisions 
with a semi-transparent EC. We set s — 3 and use a 
composite tt pulse. In order to get a superposition with 
three equal weights, we select an EC transparency of 1/3 
for the first collision (/? = 0.34) and 1/2 for the second, 
setting j3 = 0.45 after the first phase inversion. The 
sequence duration is 4.4 ms. The Wigner function of the 
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FIG. 12. Final state Wigner function W(£) for the generation 
of a three-component cat. See text for the conditions. 



FIG. 13. Final state Wigner function W(£) for the realistic 
state synthesis. See text for the conditions. 



final state is plotted in figure 12 The fidelity with respect 
to an equal weight superposition of coherent components 
centered at —0.3, 3.2 and 6.7 is 69%. 

Finally, we have simulated the state synthesis pre- 
sented in Section [V| leading to a MFSS of four coherent 
components. The pulses addressing the \h, 0) —> |g,0) 
transition can spuriously affect the \h,n) levels with 
n > 1. The frequency separation between the addressed 
and spurious transitions is quite large in this case 
for n = 1). We use thus simply an optimized square 
pulse, performing the required level mixing on the ad- 
dressed transition and a 4tt pulse on the closest spuri- 
ous transition. The maximum pulse duration involved 
in the sequence is 80 /is. For the interrogation of the 
dressed states, we use the optimized composite pulses. 
The tweezers operations are performed with a j3 — 0.6 
translation per step. The total duration of the full syn- 
thesis sequence is thus 2.9 ms. 



Figure 13 presents the Wigner function of the gener- 
ated MFSS. It is visually very similar to the ideal MFSS 
Wigner function presented in Fig. 11 (note the different 



color scales) . The fidelity with respect to the target state 
is 59%. In fact, due to the fast twezers operations, the 
final amplitudes of the coherent components differ by up 
to 0.5 from the target ones. By optimizing these ampli- 
tudes in the reference state, we get a more faithful fidelity 
of 71.5%. This value shows that a complex state synthe- 
sis operation is within reach of the planned experimental 
setup. 



VII. CONCLUSION AND PERSPECTIVES 

We have analyzed the Quantum Zeno dynamics tak- 
ing place when the photon field in a high-finesse cavity 
undergoes frequent interactions with atoms, that probe 



its state, yielding photon- number selective measurements 
or unitary kicks. A coherent classical source induces an 
evolution of the field, which remains confined in a multi- 
dimensional eigenspace of the measurement or kick. The 
quantum coherence of the evolution under the action of 
the source is preserved, the generator of the dynamics 
being the Zeno Hamiltonian, projection of the complete 
source-induced Hamiltonian onto the eigenspaces of the 
measurement or kick operators. 

The QZD evolution can be highly non-trivial. We have 
discussed in particular the generation of interesting non- 
classical states, including MFSS. We have also analyzed 
state manipulation techniques by means of phase space 
tweezers, as well as promising perspectives towards quan- 
tum state synthesis. These ideas pave the way towards 
more general phase space tailoring and 'molding' of quan- 
tum states, which will be of a great interest for the ex- 
ploration of the quantum-to-classical transition and for 
the study of non-trivial decoherence mechanisms. 

We have focused in this paper on the QZD induced by 
a classical resonant source acting on the cavity. Other 
evolution Hamiltonians could be envisioned, such as a mi- 
cromaser evolution [44] produced by fast resonant atoms 
crossing the cavity between the interrogation pulses per- 
formed on the atom at rest in the mode. The principle 
of the method could also be translated in the language 
of any spin and spring system. In particular, QZD could 
be implemented using this method in ion traps [45 or in 
circuit QED [3J] with superconducting artificial atoms. 

In conclusion, a state-of-the-art experiment appears to 
be feasible in microwave cavity QED. It would be the 
first experimental demonstration of the quantum Zeno 
dynamics. 
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